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Abstract 

We introduce a Mellin transform of functions which live on all of R and 
discuss its applications to the limiting theory of Bellman-Harris processes, 
and specifically Luria-Delbriick processes. More precisely, we calculate the 
life-time distribution of particles in a Bellman-Harris process from their 
first-generation offspring and limiting distributions, and prove a formula 
for the Laplace transform of the distribution of types in a Luria-Delbriick 
process in the Mittag-Lefner limit. As a by-product, we show how to 
easily derive the (classical) Mellin transforms of certain stable probability 
distributions from their Fourier transform. 

AMS (2000) subject classification. 44A05, 60J80. 
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1 Introduction 



In this note, we shall concern ourselves with some applications of the Mcllin 
transform to the theory of branching processes; see Sections 4 and 5 for further 
motivation. We start with a definition of the Mellin transform for functions 
that live on all of R. Next, we give an overwiew of the classical definition of 
the Mellin transform and some of its properties. This is material which, except 
perhaps for the Plancherel formula, is well-known and might as well have fitted 
into an appendix. We prove our main theorem (which is actually nothing more 
than a formula to relate the Mellin transform of a function to its Laplace trans- 
form) in Section 3. The final two sections are applications of the main theorem. 

To begin, let X and Y be two R- valued random variables such that P(X < 
x) =: T(x) and P(Y < y) =: G(y). Consider their product XY. Then, by total 
probability, 

/oo 
P(XY < z\Y G dy) P(Y G dy) 
-OO 

,0" 

= / P(X>zY- 1 \Yedy)dg(y) 

J —oo 

+ P(X< zY- x \Y G dy) dG{y) + P(0 < z) ■ g , 
Jo+ 

if Q charges {0}. It follows that if X and Y are independent, 

/0~ />oo 
(1 - Tizy- 1 )) dG(y) + / T(zy- X ) dG(y) + P(0 < z) ■ g , 
-oo JQ+ 

and if they both have a density, then 

/ © g(z) := (°° [/(zy-^giy) + f{- Z y^) g {-y)] ^ 

Jo y 

is the density of the product XY. f © g is called the Mellin convolution of / 
and g. We would like to calculate its Mellin transform (see the following section 
for the definition and a number of properties of the Mcllin transform). To that 
end, we define 

(/©<?)+(*) = /©<?(*)• 4+ (*) (i) 

and 

(f ® g)-(z) = f ® g(-z) ■ I K -(-z) = f ® g(-z) ■ I M+ (z) . (2) 

Then both (f®g) + and (f®g)~ live on R + . Now we define the Mellin transform 
of / © g as the pair 

M[f®g(x);s] := (M(f ® g)+(s),M(f ® g)-(s)) , 



3 



of Mellin transforms of / © g + and f ® g~, and similarly for any other func- 
tion that lives on R. (We reserve the symbol Aif for the Mellin transform of 
functions / that live on R + .) Now 

/>oo />oo 

M(f®g)+( S )= (f®g) + (z)z s - 1 dz= f®g(z)z s - 1 dz 
Jo Jo 

/ [f{zy- l )g{y) + f{-zy- l )g{~y)} z^ 1 dz 

o Jo V 



foo poo / \ a— 1 

' y s g(y) J o /tor 1 ) ( -J d 



z\ dy 



2// y 

s-1 



so that 



(/ © = Mf + (s) Mg + (s) + Mf-(s) Mg-(s) . 

Similarly, one finds 

M(f® g)-(s) = Mf + (s) Mg- (s) + Mf-(s) Mg + (s) . 

Define a product of two Mellin transforms according to 

M[f{xy,a]QM\g(x);s] (3) 
:= (Mf + ( S )Mg + ( S )+Mf-( S )Mg-( S ),Mf + ( S )Mg-( S )+Mf-(s)Mg + (s)) . 

We have proved 

Theorem 1 Let X and Y be random variables with densities f and g, respec- 
tively. Then 

M[f®g(x);s}=M[f(x);s]QM[g(xy,s}. (4) 

This is well-known for functions that live on R + . The multiplication is the 
same as the multiplication rule for hyperbolic numbers [S]; it remains to sec 
whether this observation has more than curiosity value. 



2 The Mellin Transform 

Let / be a complex-valued function on R + U {0} such that the integral 

P oo 

M[f(x);s] := / f(x)x s - x dx 
Jo 

exists for all complex s in the fundamental strip (a, /3) := {s € C : a < 3?(s) < 
/?}. Then M[f(x); s] is called the Mellin transform of / with respect to s. In 
what follows, we shall only be concerned with the Mellin transform of functions 
which live on R + , and consequently write JAf or Aif(s) instead of A4[f(x); s] 
when this is feasible. The Mellin transform has a number of interesting proper- 
ties. We start with 
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Lemma 1 

M[f(\x);s}=\- s M[f(x);s] (5) 



M[f(x»y,s]=»- 1 M[f(x);s/tj]. (6) 

Proof. Consider 



pOO 

M[f(\x»);s] = f(\x») X s ' 1 dx 
Jo 

poo \— s/n roo 

/ /(Ay) y s ^-V» y 1 ^- 1 dy = / f(x) x*'^ 1 dx 

Jo M Jo 



M Jo 

Now set A = 1 or \i = 1. □ 



This was easy. The following is a little harder: 

Theorem (Mellin Inversion) 2 Suppose 7 belongs to the fundamental strip 
off. Then 

f(x) = — Mf(s)x- S ds. 

2?" Jj-ioo 

We first prove the following 
Lemma 2 



Sjjf WW *-/(!)■ 



Proof of the lemma. Fix e > 0. Then, because e s is analytic in s everywhere 
in the complex plane, we have 



/ e T£ s ds = e* £ 3 ds = i 

J~{—ioo J— ioo J — c 

and so 

1 p+ioc -, ,-oo 

^- e— M/( S ) = — / e"' f(z) S- 1 dz ds 

/ e^ 3 f(e- y )e- ys dyds 

J ■y— ioo J — 00 



e-^ 2 ds = l - 
e 



2m 



7— ioo 

• 7+200 



1 pOO p^^-lOO 

= J_ / /( e -») / e^^e- 2 ^ 2 ^' 1 dsdy 
Now let e -> 0. □ 
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Proof of the theorem. By the lemma and the scaling property ([5]), we have 

1 /"7+zoo i /-7+ioo 

— / Mf(s) x~ s ds = — M[f(xy);s] ds = f(xy)\ = f(x) , 

which concludes the proof. □ 

The following theorem is in a similar spirit. 

Theorem (Plancherel Formula) 3 Suppose 7 belongs to the fundamental strip 
of both f and g. Then 



1 



7+zoo 



Mf(s) Mg(s)ds = / f(x)g(x)x 2j - 1 dx. 

7 — 'too 



Proof. As for the proof of Lemma (J2J) , one checks that 



1 

2iri 



7+ioo 



M/(s)ds = /(l) 

7 — zoo 



and 

— / Mf(s)x- S ds = — Mf(s)(~) ds = x-^f(l/x) 



Then 

1 

2?ri 



e" 8 s Mf(s)Mg(s)ds 

7— lOO 

^ /'7+^oc /*oo 



e^ 2s2 / f(x)x s ' 1 dxl y^g(l/y)y s - 1 dyds 
Jo Jo 



2ni 



2 2 

e7 re . 



7 — ZOO 



OO /'OO 



2?ri 



1 /*oo poo /«7+zoo 

— / f( e ~ U ) I e-^gJF^ I e^s'e- 2 ™ 8 ^-^ 2 ™^ 1 dsdudv 

7 — ioo 



1 />OC /• OO 

27T£ /_„, /_„ 



e- 27t ' 5 (e- t ') / /( e - 27re «-") e"™ dudv. 

Now let e — > and change variables. □ 
We now continue with the derivation of 



G 



3 The Main Theorem 

The gist of the paper is the following 

Theorem 4 Let f be a function on R with bilateral Laplace transform 



/CO 
e- ux f{x) dx 
-oo 



and Mellin transform M[f{x);s] = (Mf + (s),Mf (s)) . Define 

a+ := sup {a E R : lim e ax f+(x) = 0} (7) 

x— >oo 

and 

a~ := sup {(3 e R : lim e ax f-(x) = 0} . (8) 

Suppose that a + and a~ are strictly larger than zero, and choose for [3 any 
number in R~ such that < —f3 < min(o! + , a~). Then we have, for s > and 
in the fundamental strip of f + and f~, 



Mf + {s) = ^ v{u){-u)-°du 

2m Jp-ioo 



(9) 



and 



Furthermore, 



T(s) ff 3+ico 

M r{s) = ^7 <p(-u){-u)~ s du . (10) 

/ [3—ioG 



2ni 



¥>{u) = — / Mf+{s) T(l - s) u*- 1 ds 

-i p+ioo 

+ 55 jL. W-Mr(i -*)(-»)-' 



(11) 



/or any 7 £/iaf belongs to the fundamental strips of both f + and f . 

Proof. We first prove the theorem under the assumption that / is concentrated 
on R + , such that / = /+, and we may write Mf for its Mcllin transform. 
Choose /?:=/?+ as required. By Laplace inversion, we have 



1 rp+iao 

f(x) = — <p(u)e ux du 

2m JfJ-100 

and so 

2niMf(s)= / <p(u)e ux dux s - 1 dx. 

JO JQ-ioo 
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The standard trick here is often to first convolve / with some 
The Laplace transform of the latter is e £U / 4 , and 



/>oo pp+too 

/ / e eu2/i ip{u)e ax dux 3 - 1 dx 

JO J p-ioo 

= r +loo e^?t^r e -* x ^ dX du 

J-fi-ioo U ' S JO 

= T(s) / e eu2/4 (p(u)(~u) s du, 

J 0— too 



where interchanging the order of integration is permissible because the integrand 
becomes exponentially small if either u or x arc of the order ~ e . One can 
then undo the convolution by letting e — > 0. As for the converse, choose 7 in 
the fundamental strip of /. Then 

f( x )=2-J_. Mf(s)x- s d S , 

and 

/•GO /■'y+ioo 

/ / Mf(s)x- s dse- ux dx (12) 

JO J'y — ioo 

/•7+ioo poo /'7+ioo 

= / Mfis)^- 1 x- s e- x dxds= Mf(s)T(l-s)u s - 1 ds. 

J~{—ioo JO J~{—ioo 

To justify the interchange of the order of integration, one can Mellin convolve 
/ with (7re) -1 / 2 e^^° ex ^ / £ , whose Mellin transform is e es / 4 , and then again 
undo the convolution by letting e — > 0. This proves the theorem in case / 
is concentrated on R + . For the general case, let ip + and be the Laplace 
transforms of /+ and /~, respectively. Then 

/>oo />oo 

<p{u)= f+{x)e- ux dx+ f- (x) e ux dx = (p+ (u) + <p- (-u) , (13) 
Jo Jo 

and 



TVo\ rP+ioc p/ \ j-/}+ioo 

Mf + (s) = ^- p( u ) (_„)-« dtl _ _±i / ^-(-u)(-u)- s du, 

2m JfJ-100 lm J/3-ioo 



by what we have proved already. Now (p~ cannot have any singularities in the 
half-plane Sft(s) > — [3. We therefore can close the contour of integration in the 
second integral along a semi-circle of radius r and find, by analyticity of , 

r/3+ir fx/2 

/ tp-{-u){-u)- s du = -r ip-{re l6 ){re ie )- s d6 , 

Jd-ir J- k/2 
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which is of order 

/7r/2 pir/2 
e -rcos9 d9 < 2r 
-7T/2 



Jo 



in absolute value, and therefore tends to if s > 0. This proves ©. (fTU)) is 
proved in a similar manner, and (|11[) is an easy consequence of (fT2)) and (|13p.D 

Although we shall not use it as heavily, we note the following variant of 
Theorem [4] 

Theorem 5 Let f be a function on R with Mellin transform {Mf + ,Mf~) and 
Fourier transform 







f*(y):= e iyx f(x)dx= / e~ lyx f~ (x) dx + / e iyx f+(x) dx . (14) 

<-< — DO 

Then, for < s < 1, 

r( s ) 



Mf+(s) 



and 



Moreoi 



M (ex P (-i y ) j[°Y (y) (is) 

M/-(«) = ^»(«p(iy ) J™ r(y)y- s dy^) . (16) 



/*(!/)=< 



/■ 7+i0 ° / S7T?'\ 

y ^( S )r(l- S )exp^-j 



'7- 

/■ 7+i00 /s7rA 

^_ M/-( S )r(l- S )exp^— J y- 1 ^ (17) 



' 7 — iOO 

/or any 7 i/iai belongs to the fundamental strips of / + and f~ 
Proof. We need the following 
Lemma 3 For < s < 1 , 







e ±M, a;*- l da: = r(s)exp ( ±— ) . (18) 



dx = = —x ± 



Proof of the lemma. We have 

e ±ix e~ xu 

(i u=pi u 2 + 1 u 2 + 1 

and then, by Theorem [H 

/ cosxe Aa V 1 <£e= — — / — {-u) s du 



2« 7 7 - lco (u + A) 2 + 1 
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for arbitrary A > 0. It then follows by the Residue theorem that 



coBxe- Xx x'- 1 dx = r{s)((X-i)- s + {X+i)- s ) = T{s) 



cos ( — s arg(A — £)) 

q „ __. v „, u ,. v ,,,.,„, — 2 - i)s/2 , 

which implies 

/ cos x x s 1 dx = r(s) cos — - , 
Jo 2 

by arbitrariness of A. Similarly, one proves 

/ sinx x"~ 1 dx = T(s] sin — , 
Jo 2 

and this already implies the lemma. □ 
Proof of the theorem. Inverting the Fourier transform, we find 

/oo 
f*{y)e- i ^dy, 
-OO 

and so 

pOO pOO 

2TrMf + (s)= / f*(y)e- ixy dyx"- 1 dx 

JO J — oo 

p OO pOO />oo 

f*(-y) / e ixv X s - 1 dx dy + / f*(y) / e~ ixy X s ' 1 dx dy 
Jo Jo Jo 

\s) exp (i ^) £ f*{-y) y~ a dy + T(s) exp (-< y ) /* (y) jT* d V , 



= r( 



by Lemma[3] But this already implies the first of Equations (fl5|) - (fTT)) , because 
f*(—y) = f*(—y)- (One again checks that interchanging the order of integra- 
tion is permissible at worst after convoluting / with some (tts)" 1 / 2 e~ x2 / £ .) The 
rest are proved in a similar manner. □ 

As an application, we now calculate the Mellin transforms of those stable 
probability distributions with Fourier transform 

P(V) ■■=■■ r afi {y) ■■= exp (-bre 1 *™/ 2 ) . 

Writing e 27r6> / 2 =: £ 5 we find 

i pOO 

i>* (y) y~ s d y = / exp (-y Q c) y~ s d v 



oo p:(l-»)A 



a _1 r((l - s)/a) 



(1 - s)tt0 (1 - s)tt0 

cos 1- i sm 



2a 2a 
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because \(\ = 1. Multiplication by e tS7T / 2 and taking the real part then gives 

T(s)T((l-s)/a) (l - s) w q + sna 

Mi/j + {s) = cos — i 

an 2 a 

r(s)r((l - s)/a) . (1 - s)n(a -6) _ p + T{s)Y(l + (1 - s)/a) 



an 2a T(l + p+(l - s))T(l - p+(l - s)) ' 

if we define 

4- a — 8 , n . 
9 = ~Ja~ ' ^ 
and make use of T(s) T(l — s) = 7r(sin ns)^ 1 . p^|) is equivalent to Formula 17 
in [5] or Formula 6.8 in [BJ. Similarly, one finds 

• j I — / \ p- r( 3 )r(i + (l - a )/a) 

T(i + p (l-s)jr(i-p (i-s)j 



if we define 



P = ■ (22) 

2 a 



In particular, A / t'0 + (l) + (1) = p + + p = 1, as required. 



4 An Application to Bellman-Harris Processes 

Let {Z t }t>o be a supercritical Bellman-Harris process with offspring distribution 
{^felfeLo an d life-time distribution G. Thus, G(£) is the probability that a 
newborn individual survives at least until time t, and 717. is the probability that 
once it splits into a number Z+ of progeny, it will split into exactly k of these. 
We denote by f(s) := E(e sZ +) — J2^L nkS k the corresponding generating 
function. By supercriticality, p := /'(l) > 1, and there exists q € [0, 1) such 
that f(q) = q. We assume that G is non-lattice, that G(0 + ) = G(0) = 0, and 
that p > oo. Then there exist 'constants' Xt (the Seneta constants) such that, 
on the set of non-extinction, XtZt converges almost surely to a non-degenerate 
random variable Z whose Laplace transform ip(u) = E(e~ uZ ) satisfies 

/>oo 

^(tt)= / fo^{ue-^)dG{t), (23) 
Jo 

where /3 is the Malthusian parameter, that is, the unique root of 

/•OO 

H / e- yt dG{t) = 1 , 
Jo 

y G (0, oo). There are only a few instances where the solution of (|23p is known 
for given / and G. We turn the problem on its head: Suppose ip were (the 
Laplace transform of some function g) such that it fulfills the conditions of 
Theorem |U Then 

T(s) rP +l °° 
Mg{s) = ^-1 ^{u){-u)- s du 

Zm J 0-ioo 
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for some suitably chosen value of /3, and it follows by Equation (f25)) that 



/3+ioo poc 



Mg{s)= / /o7/>(ue-*)dG(i)(-w)- s du 



2?ri 

r(s) J/3-ioo JO 

oo />/3-Hoo 

/ / o ^(«e _t ) (-ue-*)"^-* 8 e*d(ue _t ) dG(t) 

J/3-ioo 
/3+ioc />oc 

/ o V(tt) (-u)~ s / e - *^ -1 ) dG(t) . 

I P—ioo JO 

In other words, we have for the Laplace transform of G at s, 

ip(u) (-u)-*- 1 du 

[3 — ioc 



' dG(t) 



/3+ioo 

foip(u) (-u)-"- 1 du 

[3 — too 



Example. Suppose that g(x) oc e x x K 1 for some K > and /(s) = s m . Then 

*(«)(-.)— *=^^> 

/3-ioo r ( s + !) 

e a; x dx = 2m ■ 



r(a + i)r(«) y r(* + i)r(/s) 

which eventually implies that 

r(s + k) r(m/t) 
r(s + m/c) ' 



e~ ts dG(i) 

and that G has a density: 

dG{t) _ T{mn) 



t(l _ e -f)(m-i)«-i , (24) 



rft r(mK — k) 

In case / is a polynomial f(s) — X^JLi 7r j s *' : we similarly obtain 

r(a + k) 

e~ ts dG(t)= - T M- ■ . (25) 

S ' r C*0 

It is not obvious that the function on the right-hand side of (|25p has all the 
properties of a Laplace transform (most of all, complete monotonicity) , but this 
is what we have shown. So the question is, When does the Laplace transform 
ip of the random variable Z have an exponential tail? The answer is given by 
the following 
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Theorem 6 Let F t be the PGF of particle numbers in a Bellman- Harris process 
at time t, and let f be the PGF of the corresponding first- generation offspring 
distribution. Say that f has exponential moments up to order r > if f{e u ) < oo 
for u < r, and let M t := F((l). Then F t has exponential moments up to order 
0(r/M t ). In particular, there exists r' > such that Laplace transform 

i>(u) := lim F t {e' u/Mt ) 

t — >oo 

is analytic for u > —r'. 

Proof. We can assume /(0) = 0, so Z t — > oo almost surely. Because 

{M t logF_t(e u y 



F t o exp 



V M t 



(we write F_t to denote the inverse of Ft), one readily checks that Theorem [5] 
holds true iff 

d e^ 131 M* e u 

Tu M * = F^e-^H.F.^) < £U 

(supposing u > 0) remains bounded away from zero. Now e~"'M t converges to 
something non-zero by the Kesten-Stigum theorem, so F t (e u ) — > oo for u > 
by convexity of F t . But then F t o F_t(e") = e u implies F_ t (e u ) — > 1. All that is 
left show is that e~ 0t F[ o F_ 4 (e") remains bounded. Write s instead of e u . We 
have [3] 

F t (s) = (1 - G(t))s + f f o Ft-uis) dG{u) , 
Jo 

so 

X t (s) := e-V Fl(s) (26) 

= e -/W(i-G(t))+ / / , oF u (s)e- /3u F^(s)e" /3(t - u) dG(t-u) 
Jo 

= e^(l-G(t))+|7 // °^ ) -lj ^(sJdG^t-u) 
t 



+ / X u (s)dG^(t-u) , 
Jo 

where we have introduced the measure 

Gp(t)=(i f e-? u dG(u) 



o 
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We now use the final expression in (J^HJ) to obtain 

X u {s) dG p {t - u) - / e~ Pu (1 - G(u)) dG p {t - u) 
Jo 

■ f ° Fv ^ X v (s)dG p (u-v) dG fj (t - u) 
u=o Jv=a A* 

' " " " ' ' X v (s) / dG p (u-v) dG (t - u) 



u=0 M 
t f / 

""' X v {s)dGf{t-v) 



f'oF v (s) , 



and find by induction and Fubini's theorem 



X t (s) = e- 0t (l-G(t)) + / e~P u (1 - G(u)) dU p (t - u) 

Jo 

r ( /,0 ^" ( ' S) - 1 ) X u (s)dUp(t-u), (27) 



or 

ft / f , 



*t(s) - X t (l) = jf ( /,0 ^~" (a) - l) *t-u(«) dE/^u) , (28) 
where 

oo 
i=l 

is essentially the renewal measure for G^. Suppose X t o F-t(s) does not remain 
bounded. Then we can find a sequence of values t\,t2, ■ ■ ■ tending to infinity 
such that 

/f'oF tk _ u oF_ tk (s) 



1 < lim / i Ifc ~" l^A^ _ 1 dUp(u) . 

k-Mx>J a V A 4 / 

But Zf-u/Zt — ► e - ^" almost surely [5], hence 

lim F ii _„oF_ li ( s ) = lim E (>_ t , ( S )^ ) < s^" , 
by Jensen's inequality. Therefore, 

Jo \ i" / M 

if s < e r (recall that / has exponential moments up to oreder r), because dUp(u) 
is essentially Lebesgue measure plus a term of order e~ au du for some a G (0, 0\ 
[TO] . But 1 < O(logs) is a contradiction for s sufficiently close to 1. The theo- 
rem follows. □ 
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5 Limit Laws of Luria-Delbriick Processes 



The Luria-Delbriick (LD) distribution arises as the distribution of types in a two- 
type Bellman-Harris process or, in a narrower sense, as a limiting distribution of 
types in such a process. If the life-time distribution is exponential and branching 
is binary, the theory of the LD distribution is essentially complete; see [T] . The 
same can be said if the life-time distribution of cells is exponential, but cells 
always produce a fixed number k of mutant or non- mutant progeny [5J [7J . We 
shall refer to such a process as a (1 — p, K)-Luria-Delbruck process, where p is 
the probability that upon division, a non-mutant cell produces one non-mutant 
and k mutant daughter cells. Mutants only produce k mutant progeny. The 
following theorem has been proved by Leona Schild in her diploma thesis: 

Theorem |7J 7 Denote by L n the number of non- mutants in a (1 — p, n)-Luria- 
Delbriick process when population size has reached nn + 1, and the process has 
been started from a single non-mutant individual. Then 



L (29) 



almost surely, and 



L = B 1 -^ ■ L» , (30) 

where k = n~ 1 , B is ((1 — p)x, px) -Beta distributed, and is x-biased (1 — p)- 
Mittag-Leffter. 

We explain our terms: By a (1 — p)-Mittag-Leffler distribution we mean a dis- 
tribution whose Laplace transform is the Mittag-Lefher function 

oo / \k 

^M=E r((1 ( :;; +1) - <*> 

Next, a ^-biasing of a random variable X with Laplace transform K(e~ uX ) is a 
random variable whose Laplace transform is 

E(X*e-" x ) , , 

- • (32) 



E(X") ' 

We refer to [7] for further background. Our goal here is to deduce 
Theorem |2j 8 The Laplace transform of L is 



i=0 



E (e -t) = r MS (- ) m _ p) + x y (33) 



(which has been proved in [2] in a rather indirect manner) directly from Theorem 
7. We need the following 
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Lemma 4 Suppose Y is (1 — p)-Mittag-Leffter distributed. Then its Mellin 
transform is 

[) r((l-p)(,s -i) + i)- 

Proof. First observe that 



r((i- P )(fc + i) + i) 



r((i- P )* + i) 

which implies that E\- p is analytic on all of C. Also 

dC, 



T(z) 2m J n C 

where Tt is a Hankel contour encircling the negative axis in counterclockwise 
direction. We fix e > 0, and choose for H = H(e) a vaguely lollipop-shaped 
figure as follows: It runs along (in fact, just below) the negative axis from — oo to 
— e, runs around the origin on a circle of radius s in counterclockwise direction, 
and then returns to — oo just above the negative axis. Then we have 



dt = / -r, r^—, —r, r^— dr 



d0 + mTuTTi .vn-^tj-ii dr , 



— 7T ^ 



(l—p)k+l e i6{l—p)k+i0 I r (l~p)k+l e i-rr(l-p)k+iTT 



and the three integrands are ~ e ^ p ) fc in order of magnitude. It follows that 
if | it | < £ 1_p , the sequence of functions 



j, ( 1 — p) k + 1 g — i 7T ( 1 — p) k — i 7r 
fe=0 



is uniformly convergent in r, and that of 

n 

E( 



£ (l-p)fc+lgi6»(l-p)fc+ze 
fc=0 



is uniformly convergent in 0. It is therefore permissible to interchange the order 
of summation and integration in the following chain of equations, and we find 

C 



1 _ _ c c 



1 f e < " {-u) k 1 f e c 

^Jn~hc^ dC ^^JnCT^ dC - (34) 



Now the right-hand side of Equation (|34[) defines an analytic function in u except 
possibly for those it which belong to the zero set TCo(e) := {u E C : C + U C P = 
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for some £ G H(e)} of the denominator. But one readily checks that 7io(s) is 
a Hankel contour around the positive axis that winds around the origin along a 
segment of a circle of radius £ 1_p , so that if we choose some e > e, we have 



H(e) C + < P Jn{e) C + < 

on a (connected) segment of the open torus (ugC: e 1_p < \u\ < e 1 ~ p }. By 
analytic continuation, then, (|34l) holds for all u € C \ TLo(e), and in particular 
holds for all u for which 5R(w) < — e 1 ~ p . By Theorem[3J therefore, we can write 

if j3 < — e 1_p . Now one may interchange the order of integration (if necessary 
again after convoluting with y^rree~ x / e ) and then apply the Residue theorem 
to find 

T(s) f r P+lcc e c 
MY{s) = -^± / -——(-u)- s duaK 



If e c 



T(s) 



((l-p)( s -l) + l) ' 

which again follows from Hankel's representation of the reciprocal of the Gamma 
function. The proof of the lemma is complete. □ 

Next we need 

Lemma 5 Suppose B is ((1 — p)x, pxj-Beta distributed. Then the Mellin 
Transform of B l ~ p is 

MB^(s) = T{X) r((l -,)(«-! -fx)) 

[) r((i-p)x)r((i- P )(*-i) + x)- 

Proof. We have 

P(B <x) = T [ X \ [ X y^-e)*-\l - y)^ dy . 

T((l - p)x)T{px) Jo 



Then P(B 1 - p < x) = P(B < x^-^ '), and 

y (i-p)*-i(l-y)p*-idyx'- 1 da; 



1 d ^-" r 



o dx J 



^ / ,/' ">* '(1 ,/)'—'. /,/,/' '"<" l '</„ 



, r((i-p)(«-i) + x) 
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which already finishes the proof. □ 

It follows from Lemma 2] that the Mellin transform of a x-biased Mittag- 
Leffler distribution is 

Together with Lemma [5] and Theorem [1] this implies that the Mellin transform 
of L satisfies 

r(s + x- 1) 

= r«i -V i>' + ») • 

Specialising now to s = i + 1 for arbitrary integer i > 0, one has for the Laplace 
transform of L, 

E(e' uL ) ^2, r(z + x) (-u) 4 



E 



r(x) ^r(x)r((i- P )i + x) i 



oo % 



V TT - + - ' ~ 1 (~ M ) = y 



J r((i-p)« + x) ^ V * yr(z(i- P ) + x) ' 

This was to be proved. □ 
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